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ABSTRACT 

String theory abounds with hght scalar fields (the dilaton and various moduli) which 
create a host of observational problems, and notably some serious cosmological difficulties 
similar to the ones associated with the Polonyi field in the earliest versions of spontaneously 
broken supergravity. We show that all these problems are naturally avoided if a recently 
introduced mechanism [16] for fixing the vacuum expectation values of the dilaton and/or 
moduli is at work. We study both the classical evolution and the quantum fiuctuations of 
such scalar fields during a primordial infiationary era and find that the results are naturally 
compatible with observational facts. In this model, dilatons or moduli within a very wide 
range of masses (which includes the SUSY-breaking favored ~ 1 TeV value and extends up to 
the Planck scale) qualify to define a novel type of essentially stable ultra-weakly interacting 
massive particles able to provide enough mass density to close the universe. 



I. INTRODUCTION 

Superstring unification [1] and inflationary scenario [2] have been, arguably, the most 
influential ideas in particle physics and cosmology over the last decade. Attempts at combin- 
ing these ideas and constructing superstring-based inflationary models have introduced some 
interesting new ideas [3], but have, however, encountered serious difficulties [4-7]. A speciflc 
source of difficulties is the existence of massless scalar flelds (the moduli) having only gravita- 
tional strength couplings to ordinary matter. Among the moduli flelds, the dilaton $, which 
invariably accompanies the graviton in all superstring models, plays a special role. In the 
presence of a (tree-level coupled) dilaton, the cosmological evolution is drastically different 
from that in Einstein's gravity. In particular, instead of driving an exponential inflationary 
expansion, a constant vacuum energy drives the dilaton towards large negative values (corre- 
sponding to weak couplings) , while the universe expands only as a small power of time. Even 
apart from inflation, a massless dilaton, or moduli, held (for simplicity, as these flelds share 
many characteristics, we shall refer to any of them as " the dilaton" and denote it as $) can 
cause a host of other cosmological problems. During the most recent, matter-dominated epoch 
of the universe, such held is necessarily time-dependent. The masses of elementary particles 
and their couplings all depend on and the cosmological time- variation of the dilaton runs 
into a violent conflict with observations, not to mention unacceptably large violations of the 
equivalence principle. 

These problems are usually addressed by assuming that the dilaton develops a poten- 
tial, so that one recovers the standard Einstein's gravity after $ settles at the minimum of 
the potential. The potential could originate from non-perturbative effects, such as gaugino 
condensation. The existence of a potential for $ entails, however, new difficulties. On the 
one hand, it has been argued [7] that the minimum of the nonperturbative potential is too 
shallow to conflne the dilaton without flne-tuning the initial conditions (see however [8]). On 
the other hand, any potential V{^) for a weakly coupled fleld resurrects the Polonyi problem 
[9-12]: either the energy stored in the coherent oscillations of the vacuum expectation value 
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(VEV) of ^» does not dissipate before now and exceeds the critical density needed to cfose 
the universe, or the field decays before now , and thereby generically produces an excessive 
amount of entropy. In the case of the nonperturbative potentials suggested by current SUSY 
breaking models in string theory, the slow decay rate of the moduli fields leaves the universe 
in a radiation-dominated era at a temperature which is generically (i.e. under naturalness 
assumptions for the couplings) much too low to be consistent with primordial nucleosynthesis 
[13-15]. 

In the present paper, we show that all the above difiaculties associated with the dilaton 
(and moduli) fields are naturally avoided if the new mechanism introduced in Ref. [16] for 
fixing $ is at work [17]. The main idea of [16] was to exploit the fact that string loop 
effects, associated with worldsheets of arbitrary genus in intermediate string states, naturally 
generate some non-monotonic dependence upon $ of the various couplings of $ to the other 
fields. Under the assumption that the different coupling functions -Ba($) have extrema at some 
common value of $ = $0) it has been shown that interactions with massive particles in an 
expanding universe drive the dilaton towards the value $o ^tt which it decouples from matter. 
All deviations from general relativity in this model have been estimated to be extremely 
small at the present cosmological epoch thereby naturally reconciling a massless dilaton with 
existing observational data. The assumption of coincident maxima can be satisfied in a 
technically natural way if there exists a discrete symmetry, e.g., $ — > — $. Such a symmetry 
would guarantee that all couplings have extrema at = (the only additional necessary 
assumption being that these extrema lead to minima rather than maxima of the masses as 
functions of $). We note that precisely the discrete symmetry $ — > — $ (or — > l/e*^) is 
known to hold for some of the moduli fields (T-duality), and has been conjectured to hold for 
the dilaton proper (S-duality: Qs — > l/fi's, where Qs = e* is the string coupling). 

Here we extend the analysis of Ref. [16] to inflationary models. It will be shown that 
inflation is extremely efficient in driving a homogeneous fleld $ to $o (Section II). At the 
same time, inflation is known to generate signiflcant quantum fluctuations in gravitational 
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and other fields with a very wide range of wavelengths. Fluctuations of the dilaton on co- 
moving scales smaller than the present horizon are potentially dangerous because they are 
not damped by the mechanism of Ref. [16]. One might worry that quantum fluctuations 
could reintroduce the Polonyi-moduli problem in a different form. However, in Sec. Ill we 
shall calculate the fluctuation spectrum and show that the predicted dilaton fluctuations are 
well below the observational constraints. We next consider, in Sec. IV, the possibility that, in 
addition to non-trivial matter-coupling functions with extrema at $ = ^»0) the dilaton also has 
a potential with a minimum at $o- (This may be enforced by the same discrete symmetry). In 
this case the dilatons are massive, and can potentially run into conflict with observations, e.g. 
by overdosing the universe or by generating an excessive flux of 7-rays due to dilaton decays. 
We shall see, however, that in our model the constraints on the dilaton mass derived in [11-13], 
[18] can be substantially relaxed, due to the very weak couplings of the dilaton. We flnd that 
dilatons with a very large range of masses (which includes the suggested SUSY-breaking value 
~ TeV and extends up to the Planck scale) can qualify to deflne a new type of (essentially) 
stable WIMP able to provide enough mass density to close the universe. 

II. THE EFFECTIVE ACTION 

Up till now, most of the analysis of superstring cosmology has been based on the tree-level 
effective action, corresponding to the lowest order in the string loop expansion, 

Stree = j dt^X^ge'^'' {{a'y^ R + 4V2$ - A{y^f] + £^atter}. (2.1) 

Here, $ is the dilaton, and the matter Lagrangian includes fermions, gauge and Higgs flelds, 
and in particular the 'inflaton' scalar fleld % whose potential V{x) drives the inflation, 

>C„.atter = " ^^V' - ^(Vx)' - V{x) + •••• (2-2) 

Hats in Eqs.(2.1), (2.2) indicate that the corresponding fields are taken in the 'string frame', 
that is, in the cr-model formulation of string theory. 
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The string coupling Qg, which plays the role of the expansion parameter in the string 
loop expansion, is determined by the expectation value of the dilaton, Qs = e'^ . The tree- level 
action (2.1) is proportional to g~^, resulting in a universal, multiplicative coupling of the 
dilaton to all other fields. With higher orders in the loop expansion taken into account, we 
expect the common factor e~^* to be replaced by several coupling functions Ba(^) multiplying 
different terms in (2.1), (2.2). In particular, the effective action for the graviton-dilaton- 
inflaton sector is expected to be of the form 

S = jd'x^g\^^R^^[AVH-A{V^f]-\B^m^^^^ (2.3) 

where the functions Ba{^) admit a series expansion 

B,{^) = + + c^^^e^^ + (2.4) 

and a similar expansion for F(x, $). In the case of the moduli fields (by contrast with the 
four-dimensional dilaton proper) the effective action has also the generic form (2.3), the only 
difference being that the non-trivial ^-dependence is absent at tree level, and arises at one 
loop and beyond. 

A more convenient form of the action can be obtained by a conformal transformation 
from the 'string-frame' metric g^i^ to the 'Einstein-frame' metric 

g^, = CBg{^)g^,, (2.5) 

and by replacing the dilaton field $ by the variable 



1/2 



This gives 

S = I d'^x^ 



(2.6) 



(2.7) 



where we have defined 

X = C-'/'x, (2.8a) 



F{v) = B^i^)/B,i^), (2.86) 
V{x, <p) = C-^B;\^)V{x, (2.8c) 

The constant C in Eq.(2.5) is chosen so that the string units coincide with Einstein units 
at the present cosmological epoch, CBg{^Q) = 1, and the constant q in (2.7) is related to 
the bare gravitational constant G, q — AttG = Ca'/A. As shown in [16], G is numerically 
nearly identical to the observed Newtonian gravitational constant so that q — An /nip with 
rup = 1.22 X 10^9 GeV. 

The minimal condition required for the mechanism of Ref. [16] to work is that all -Ba(^) 
should have an extremum at the same $ = $o- When formulated within the context of 
inflationary models and in terms of the rescaled fields (p and X) this leads to requiring that 
the potential V(x, ip) in Eq.(2.7) has a minimum (as a function of (p) at ip = po for any fixed 
value of X- Here, (po = <^($o)- A simple toy model which satisfies this condition is the case 
where higher loops are supposed to respect the tree-level universality of the dilaton couplings. 
In this case all functions Ba{^) in (2.3) are identical, 

S„($)=S($), (2.9a) 

and 

V{X,^) = B{^)V{X). (2.96) 
Then, in the Einstein-frame action (2.7), F{(p) = 1 and the potential has a factorized form. 



Tq^ - ^(^^)' - ^(^X)' - B-\p)V{x) 



(2.10) 



In this model ip is attracted (both during inflation and the subsequent matter-dominated 
era(s) ) toward the maxima of B{(p). Although the strong universality condition (2.9) may 
be too restrictive, it is useful to have in mind the action (2.10) which provides a very simple 
model containing most of the essential physics of the more general model (2.7). 
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III. CLASSICAL EVOLUTION 

We shall assume, as it is usually done in inflationary models, that the potential V{x, 
in Eq.(2.7) has a minimum with ^(xo, ^^o) ~ and that it has a 'slow-roll' region where it is 
a slowly varying function of x, 

V^'^ < 12qV^, (3.1a) 

where we recall that q = An/m^. At the same time, the coupling functions Ba and the 
potential V are not expected to be slowly-varying functions of (p. For example, in the strong 
universality model (2.9) we expect B{(Pq) ~ 1 and 

« = -B"{<po)/B{<po) ~ 1 (3.2) 

at the maximum of B{(p). Note that, contrary to x which has the usual dimension of mass, 
(fi is a, dimensionless variable whose expected range of variation is of order unity. 

For the initial conditions of the universe, we shall assume that the fields (f and x ^ire 
displaced from the minimum ofV{x, <fi) at {xo; fo} but are in its basin of attraction, at least 
in some parts of the universe [19]. Then it is easy to see, qualitatively, how the cosmo logical 
evolution will proceed. The field <f, which corresponds to the steep direction in V{x,<f), 
will evolve on a much faster time scale than x- It will start oscillating about (p = (po, 
these oscillations will be damped by the expansion of the universe, and the universe will 
quickly settle into a quasi-exponential infiation driven by the potential energy of the field 
X) ^(x) = ^iXjfo)- The damping of (^-oscillations during inflation is very efficient, and we 
expect that by the end of inflation the dilaton will be very close to ipo. 

To describe this quantitatively, let us consider the field equations for x and (p 

V(F(^)Vx)-|^nX,^)=0, (3.3) 
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During the slow-roll phase of inflation, the universe can be locally described by a flat 
Robertson- Walker metric, 

ds'^ = -dt"^ + a'^{t)dx^, (3.5) 
with the expansion rate H = d/a given by 

H^ = l [2gy (x, ^) + 0V qF{v)x'] « ^qV{x, ^o)- (3.6) 

The spacetime is approximately de Sitter, with the curvature 

R ^ 12H^ ^ 8qV(x, ipo). (3.7) 

To study the approach of (p to (po, we expand Vlxjf) in powers of {(p — po). For 
notational consistency with [16], it is convenient to denote by Pi the (positive) dimensionless 
parameter measuring the curvature (with respect to (p) of the inflationary mass scale around 
the minimum (pg: 

A,(x, ^) ^ V^l\x. ^) ^ A,(x, ^o)[l + \^{^ - ^^f\. (3.8) 

Here, the index i stands for 'inflation'. In the simple model (2.9), in which the potential is 
factorized, y(x, <^) = B~^{(p))'V{x)-, one has /?j = k/4, where k was introduced in equation 
(3.2) above. In the general case, the value of /?i depends on the physics determining the mass 
scale Aj. As discussed in Ref. [16] (see Eq. (4.6b) there), if the hierarchy A^ <^ (which 
seems necessary in inflationary models) is due to nonperturbative effects, one expects to have 
A ~ ln(Astring/Ai)>l . 

Substituting this in Eq.(3.4) and using (3.6), (3.7), we obtain 

(V^ - iiR)b^ = 0, (3.9) 



where S(p = (p — (po and 
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(3.10) 



We note that Eq.(3.9) has exactly the form of that for a massless, non-minimally coupled 
field. The curvature is a slowly-varying function of time during infiation, and the effect 
of the curvature term in (3.9) is similar to that of a positive mass squared term rri^ ~ 
H'^ for the fluctuations of ip. Inflation is followed by a x-dominated expansion when the 
inflaton field x oscillates about the minimum of its potential (say V ~ ^m'^{(p)x^)- During 
this x-dominated period, S(f will approximately satisfy, after taking an average over the x 
oscillations, an equation of the same form as (3.9) but with a different value of ^-^ = 
(here (3^ = [d'^\nm^{(p) /d(p'^]ip=^Q measures the (/7-curvature of the mass of the x field; see 
the Appendix for a more exact equation satisfied by Sip). 

Neglecting the spatial gradients of (p (which are rapidly suppressed by the cosmological 
expansion) we can finally rewrite (3.9) in the form 

Sip + SHSip + U^iH^Sip = 0. (3.11) 

During the slow roll period, the field (p{t) changes at a much faster rate than H{t), and 
we can solve (3.11) using a WKB-type ansatz, 

Sip^e^^'\ (3.12) 

and assuming \W\ <^ W^. Substitution of (3.12) into (3.11) gives a quadratic equation for W 
with the solution 

W{t) = ± yf^j H{t). (3.13) 

Alternatively we can follow the method of Ref. [20] and write the equation describing the 
evolution of (p in terms of the parameter p measuring the number of e-foldings: 

p = ln^= Tm, (3.14) 

where U denotes the time at the onset of inflation. The latter equation can be written in 
all eras of interest, inflation, x-dominated (in the approximation where one averages over 
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the X oscillations), radiation-dominated, and matter dominated. When neglecting (<^p)^ with 
respect to 1, this equation reads 

'^-^l^ + (1 - A)(^; = -(1 - 3A)/?(^ - ^o) (3.15) 

where A is the ratio between pressure and energy density (i.e. —1, or 1/3 in vacuum-, or 
matter-, and radiation-dominated eras, respectively), and where j3 is the parameter measuring 
the curvature of the (/^-dependent mass scale driving the evolution of in the corresponding 
era (e.g. (3i in the inflationary era and (5^ in the x-oscillation dominated era). 

The equation (3.15) is that of a damped harmonic oscillator. The critical value of (3 
separating the overdamped-type solution from the damped-oscillation-type one is jSc = 3/8 
(in both the vacuum-dominated and the matter-dominated cases). During inflation, the 
approach of (p toward <^o is oscillatory if Pi > /?c, i.e. < 0, where we define 

= 6(/?e - A) = ^ - 12^^ (3.16) 

From what we said above, we expect to be in this regime (v = ilvl). Then 

S(p = Ae-^P/^cos(\u\p + S), (3.17) 

where A and S are constants. The number of e-foldings during inflation is Pi ^ 65, so that if 
initially \Sip\ ~ 1, i.e. A 1, then, by the end of inflation, 

\5ip\<e-^P^/^<10-^^. (3.18) 

During x-domination the average pressure of the oscillating x-fleld is nearly zero, the 
expansion law is a{t) oc t^/^, and the approach of (p toward (po is described by Eq. (3.15) with 
A = and f3 = f3^. This yields, when (3^ > an additional factor of attraction of 6(p> toward 
zero of order e"^^'^/^, where is the number of e-foldings during x-domination. [Note that 
in the simple model (2.10) one has = |k = 2/?j.] The energy of x-oscillations eventually 
thermalizes, and the universe enters the radiation era with a{t) oc t^/^. Since A = 1/3 or 
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R = then, we see either from (3.15) or from (3.9) that 6(p essentially stops evolving during 
the radiation era independently of the value of (3. During the subsequent matter-dominated 
era, 5(p will be further attracted toward (po by an additional factor e"^^'"/'^, where Pm ~ 9 is 
the number of e-foldings separating us from the end of the radiation era. Finally, the present 
value of 6ip is expected to be 

\5if\< 10-^^. (3.19) 

Note that the estimate (3.19) is independent of the precise values of (3i and as long as they 
are both > 3/8. A nonzero value of 5(fi causes a number of potentially observable deviations 
from general relativity [16]. However, all observable non-Einsteinian effects are proportional 
to the square of Sip. The present observational bounds (from equivalence principle tests) are, 
within the context of the model (2.9), 

^I'^V'lobs^Sx 10-^ (3.20) 

Even if we were to assume that Pi, or equivalently ^i, is anomalously smaller than unity, it 
would only be in the extreme case where < 10~^ that the attraction factor due to inflation 
~ e~^^*P* would not be much smaller than unity. We see that inflation is extremely efficient 
in driving a homogeneous, classical fleld (p to (fQ. In the case , considered below, where the 
held (p has a potential (sharing the discrete symmetry of the coupling functions Ba{(p) ), we 
conclude that we have here a natural , non flne-tuned, solution to the Polonyi-moduli problem 
as the VEV of (p is left, at the end of inflation, very precisely at the place where it stores 
no potential energy. However, as the change in the equation of state at the end of inflation 
can result in copious creation of dilatons we must consider whether this can regenerate a non 
trivial quasi-classical VEV for (p. This will be discussed in the next section. 

IV. QUANTUM CREATION OF DILATONS 

Particle creation during and shortly after inflation can be studied using the standard 
methods of quantum field theory in curved spacetime [21]. For a massless scalar field with 



coupling to the curvature, as in Eq.(3.9), this has been done by Ford [22] in the hmit where 
the coupling to the curvature is nearly conformal, — 1/6| ^ 1. He assumed also that the 
X-dominated period is very short, so that inflation is followed by thermalization in about one 
Hubble time. He found that the energy density of created particles at the end of inflation 

{t = t*) is 

p^{U)^10-\6Ci-lfHt (4.1) 

where is the expansion rate at t = t*. 

We have calculated the energy density and the spectrum of dilatons without assuming 
{^i — 1/6) to be small and without assuming that the x-dominated period is necessarily short. 
We found that the matching to a x-dominated expansion brings several qualitatively new 
features but does not change drastically the overall quantitative results of a matching to a 
radiation-dominated era. To keep things simple we discuss in the text only the latter case. 
The details of our calculation are given in the Appendix (which contains also a brief discussion 
of the matching to a x-dominated expansion), and here we shall only state the results. 

The spectrum is expressed in terms of the co-moving wave number k, which is equal to the 
physical momentum of the wave at t — (we set a(t*) = 1)- At later times the momentum 
is p{t) = k/a{t), and the wavelength is X{t) = {27r/k)a{t). We find that the energy spectrum 
of dilatons is peaked at A; ~ i?* and is exponentially suppressed for k ^ H^. In the limit of 
long wavelengths, k -C i?*, 

for u real (i.e. < 3/16 in Eq. (3.16)) and 

for v = ilv], where in the latter case we assumed for simplicity that exp(— 7r|i^|) -C 1. [Note 
that = 1/6 corresponds to = 1/2 for which Eq. (4.2) indeed predicts no particle produc- 
tion.] 



When inflation is foUowed by x-domination, the dilaton spectrum contains a second peak 
at /c ~ which contributes roughly the same number density as the peak at /c ~ and 
an energy density differing by a factor ~ m^/H^. By integrating the spectra (4.2) or (4.3) 
up to /cmax ~ wc find (roughly independently of v as long as l^- — 1/2|>1) that the total 
dilaton energy density is of order 



p^{t) ^ lQ-'Hta-\t). (4.4) 

In other words, we find that when 6^^ — 1 becomes > 1, the factor (6^j — 1)^ present in the 
nearly conformal case (4.1) saturates to something of order unity. 

Using the approximate conservation of the comoving entropy oc NT^a^ and Hubble's law 
at thermalization [H^ ~ A/"* ^^T^/m^) we can write for the present energy density in massless 
dilatons 

Here pc is the critical density, H is the number of spin degrees of freedom in radiation and 
relativistic particles (at present), H* is its value at (i.e. at thermalization), VLj. = pr/pc = 
4 X lO^^h^^, pr is the radiation density, /i ~ 1/2 is the Hubble parameter, and nip is the 
Planck mass. The characteristic wavelength of dilatons for the peak at /c ~ ff* (when > 1/12 
[23]) is 

Ac - 27tH-^Z^ - 4(mp/if*)^/Vm, (4.6) 

where ~ Affi Af-^H^rUp)^/^ /T, with T = 2.74i^ = (0.83mm)-\ is the redshift at t = U. 
The presence of a x-oscillatory era would add a (possibly overlapping) second peak with 
characteristic wavelength differing by a factor ~ H^/m^. 

Waves of wavelength smaller than the Hubble length, \{t) < t, can be treated classically, 
and their amplitude can be estimated from 



(fikit) ~ a{t) 



q dp^{t) 
k dk 



nl/2 



(4.7) 



and at the present time we find, using Eqs.(4.2), (4.3) 



IV > 0, 



(4.8) 



(4.9) 



For k ^ H^, corresponding to wavelengths ~ Ac, both of these equations give 



-34 



(4.10) 



where we have used the bound [2] < 10 ^nip on the rate of inflation. (Larger values of 
result in an excessive amount of relic gravitational waves). In the case of z/^ < 1/4 (^j > 1/6) 



longer wavelengths. This means that for the wavelengths and time scales of relevance to 
laboratory experiments the quantum-regenerated dip is many orders of magnitude below the 
observational bounds (3.20). 

A different behavior is obtained for u > 1/2 {^i < 1/6), when (pk grows towards longer 
wavelengths. The largest growth occurs for nearly minimal coupling, f« 0, z/ ~ 3/2, when 
(pk oc A;~^+^^% and on the present Hubble scale k ~ Ha ~ ffZ* with ~ 10^^ cm 



S^m.. ~ nl/\H/Tf^^{Hjmpf-^^^ 10-^10-10^^^ (10^i^,/mp)l-2^^ (4.11) 



Even for anomalously small curvature couplings ^j^lO~^, and the maximal allowed value 
of H^, the dilaton amplitude is smaller than Sip ~ 10~^, and therefore smaller than the 
observational limits (3.20) which become less stringent when k — is itself small, i.e. when 
the interaction of dilatons with matter is suppressed due to a small value of ^j. 

It should be noted that Eqs.(4.7)-(4.9) cannot be used for wavelengths longer than the 
Hubble length, where real particles cannot yet be distinguished from vacuum polarization 
effects. In particular, we cannot conclude from (4.8) that (p>k can become arbitrarily large in 
the limit of long wavelengths. To estimate the dispersion of the dilaton field on super-horizon 
scales, we calculated the quantum expectation value (v^^). The calculation is outlined in the 



[including < 0, i.e. > 3/16] the dilaton amplitude (pk oc k^ 



2 decreases towards 



Appendix, and the result is that for > 1/2 we have {<p^)it) ~ Vkffi where kn corresponds to 
the Hubble scale at time t. This indicates that super-horizon wavelengths do not significantly 
contribute to S<f. 

V. MASSIVE DILATONS 

Up till now we have been considering the case where the dilaton remains exactly massless 
at low energy. However, as we remarked above, under the assumption of a discrete (^-symmetry 
(or some other universality feature such as the one built in the model (2.9)) the existence of 
a mass term for (p does not create the usual Polonyi problem because, after inflation, the 
VEV of (fi is left very precisely pinned at the minimum of its potential. We must, however, 
investigate what constraints on the dilaton mass m^p are obtained by requiring that the present 
mass density of the quantum- generated dilatons does not exceed the critical density pc- To 
simplify the discussion, we shall consider only the case of z^^ < 1/4, when long- wavelength 
modes with k <S are unimportant, and the particle interpretation of the field S(f becomes 
valid shortly after t^. The mass and the dilaton number density should then satisfy 
the condition 

n^ = ^^^<i. (5.1) 

Pc 

To estimate n^, we note that the ratio 

r = n^/rir, (5.2) 

where rir is the density of particles with masses smaller than the temperature, remains approx- 
imately constant in the course of cosmological evolution (assuming that the dilaton lifetime 
exceeds the age of the universe, see below). At the end of infiation, we find from integrating 
the number density spectrum (see Appendix) 

n^{U) ~ l^-'^Hl. (5.3) 

Using the other equations 

n^(t*)~ACT3, (5.4) 
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(5.5) 



where is the thermahzation temperature, we get 



(5.6) 



Inserting this in (5.1) we obtain 



% ~ lO-^A/; ~'{Hjmpf/'^{m^/T)VL. 



(5.7a) 



where, as above, VLr — Pr/Pc = 4 x 10 



Numerically this reads 



3 




lOGeV 



(5.76) 



When the equality fi.^ = 1 is satisfied, dilatons dominate the mass density of the universe. 
This can happen for the whole range of masses m<^ ^ lOGeV. Let us note in particular that the 
value ~ ITeV suggested by many current SUSY breaking models [13], [14] is allowed and 
corresponds to ~ 10~^mp, i.e. ~ lO-'^^GeV. We remark also that the value ~ rup 
[24] corresponds to a "weak scale inflation" ~ lOOGeV [17] (i.e. to an intermediate scale 
reheating temperature T* ~ 3 x 10^°GeV ~ {mw'nipY^'^). Other authors have suggested the 
possibility that dilatons may provide the dark matter of the universe [18], [25]. An important 
difference of our model is that our allowed mass range is ^ lOGeV (which contains notably 
m^p ~ ITeV) and still corresponds to essentially stable dilatons, with a decay time much larger 
than the age of the universe (as is discussed next). [Usually [18], stable dilatons exist only 
for < lOOMeV because of Eq. (5.9) below.] 

Let us finally examine whether additional constraints on the dilaton mass follow from an 
eventual fliix of 7-rays resulting from the dilaton decay, (/? — > 77. This decay is described by 
the term 



Ant OC ((/? - >^o)F^^ 



(5.8) 



in the effective Lagrangian. When expressed in terms of a canonically normahzed scalar field 
fcan ~ V'/x/?: Eq. (5.8) contains a coupling constant oc l/nip. By dimensional analysis, if 
the dimensionless coefficient in front of (5.8) is of order unity, the corresponding lifetime is 

T<p ~ mj/mj. (5.9) 

In our case, however, all coupling functions are expected to have extrema ai (p = (po, and 
thus the dimensionless coefficient in front of (5.8) is oc 5(p, i.e. is exceedingly small. Then 
7-rays can only be produced in binary collisions (fip — > 77. The corresponding interaction 
term is 

AntOc((^-(^0)'i^^l, (5.10) 

with a dimensionless coefficient of order unity. The corresponding annihilation cross-section 
is extremely small [26] , 

a r-. ml/vml, (5.11) 

where v is the average velocity of dilatons. The lifetime can be found from riypavTyp ~ 1, 
which gives 

^ rrip/n^ml,. (5-12) 
The annihilation rate per unit spacetime volume is 

^~^^4-'o-, (5.13) 

where to is the present age of the universe. Hence, there is no more than one annihilation in 
the entire visible universe during its whole lifetime! 

Note that as all coupling functions have extrema at <^ = <^0) the type of dark matter our 
mechanism leads to has only exceedingly weak interactions with ordinary matter. Basically, 
its presence can be felt only through the gravitational effect of its mass. This leaves little 
hope of detecting it in laboratory experiments. For instance, for a macroscopically sizable 
mass ~ rup ~ 2 x 10~^g, the average present cosmological density of dilatons would be 
at most ~ 10~^^cm~^. 
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APPENDIX 

In this Appendix we study the spectrum of created dilatons in simple models in which 
de Sitter inflation is followed either by a radiation-dominated expansion or by a x-dominated 
one. We consider first the transition to a radiation-dominated era. Using the conformal time 
coordinate drj = dt/a{t), the corresponding metric can be written as 

ds^ = a^{r]){-dr]^ + dx'), (Al) 

a{r]) = -{H^r])-^, rj < ry*, (A2a) 

a{r]) = H^r] -f]), r] > r]^. {A.2b) 

Here, = const is the expansion rate during inflation, 77* < is the thermalization time 
when inflation ends, and f] = r]* + {H^r]^)~^. It will be convenient to set r/* = —H~^, so that 
a(?7*) = 1. The calculation in this Appendix follows the standard techniques reviewed in [14]. 

The fleld operator 5(f{x) = (p{x) — (po satisfles a massless, non-minimally coupled fleld 
equation (3.9), 

and can be expanded in terms of creation and annihilation operators, 

1/2 r j3jL 

^^^"^ = !m / [«kV.(r/)e^''- + h.c.] . (AA) 

Here, hats indicate operator quantities and should not be confused with the notation of Section 
II where they indicate quantities in the 'string frame'. The mode functions i/^kiv) satisfy the 
normalization condition 

V^fcVfc - i^ki^k = -h (-4.5) 



corresponding to [a^, aj^.,] = 5(k — k'). We shall assume that the quantum state of the dilaton 
field during the inflationary period < ?]* is the do Sitter-invariant Bunch-Davis vacuum (i.e. 
that tpkif]) ~ e~*^''/-\/2^ when -q — > — oo). The corresponding mode functions are 

Mv) = A{-riy/^Hi'\-kri), (A6) 

where 

A = (7r/4)^/2gi7r/4gi7.^/2^ 

are Hankel functions, and u is given by Eq.(3.16). 
The mode functions for the radiation-dominated period r] > rj^ are 

and Eq.(A.5) gives a normalization condition for ak and Pk, 

The coefficients ak and Pk can be determined by matching the mode functions (A. 6) and 
(A.S) and their derivatives at = ry*. The dilaton spectrum (in number density and energy 
density) can then be found from 

dn^ = 27r2a3(t) l^^l^^^^^' {A.lOa) 

The instantaneous thermalization model (A. 2) is not adequate for modes with wave- 
lengths shorter than the de Sitter horizon at r]^, k>H^. In a more realistic model, the 
transition from the vacuum to the radiation equation of state takes at least a Hubble time, 
and particle creation in such modes is exponentially suppressed, (3{k ^ H^) a; 0. In the 
opposite limit, k <^ H^, we can find and Pk using the asymptotic form of the Hankel 
functions for small values of the argument. 



For u > 0, the first term in (A. 11) is negfigible, and we find CKfc ^ —(3k and 



In deriving (A. 12) we have used the identity 

V{z)V{l-z) = -^ (A13) 

In the case of imaginary v = we shall assume for simplicity that ex.-p{—TT\v\) <^ 1. 
Then again the first term in (A. 11) can be neglected, and we obtain ~ —Pk with 

Combining (A. 12) and (A. 14) with (A. 10) we obtain Eqs.(4.2),(4.3). Let us note in passing 
that the spectrum (A. 14) is formally of the Rayleigh-Jeans form: |/3a;|rj — Trj/^, that is the 
limit of (e'^/^i^J — 1)""*^ when k -C Trj. (Similar spectra also come out of the superinflationary 
scenario of Refs. [3], [18]). By contrast, we shall see below that when inflation is followed by 
X-domination one does not get a Rayleigh- Jeans-type spectrum. 

The dispersion of the dilaton field on super- horizon scales , /c <C ry"^, can be estimated 
by calculating the expectation value 

The contribution of long wavelengths {k <^ H^) to this integral can be found using 
Eqs.(A.8),(A.12) and (A. 14) with ak ~ -Pk- This gives 

^^^^^'^ " An^H^lv-f,r J d'kk-'\Pkfsm'[k{v-v)i (A16) 
and for > we find 

It is easily seen that for u < 1/2 the integral is dominated by the upper limit (A; ~ i?*), while 
for 1/2 < ly < 3/2 the dominant contribution is given by /c ~ ?7~^, that is, by the modes of 
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wavelength comparable to the Hubble length. In the latter case the upper limit of integration 
can be extended to infinity, and we obtain a somewhat unwieldy expression 

((5(^)2) ~ -27r-3sin(i/7r)r(l - 2iy)T\iy){iy - l/2)\H^[4a{ri)]'^''-^, {A.18) 

where a{r]) is given by (A. 2b). Disregarding numerical factors, this gives 

mr)^qH'Aa{v)?''-'. {A.19) 

This is to be compared with the dilaton amplitude (pk which can be found from 
Eqs.(4.7),(4.2), 

2,^, qk^ fH,^^-'+' 



It is easily seen that ~ ip^ for k ~ r]~^. 

Let us finally briefly indicate the peculiarities of the more general case where the in- 
flationary evolution (A. 2a) is followed by a x-dominated expansion with an averaged scale 
factor (we assume <^ H{r]) <^ away from r] — r]^ and average over the oscillations at 
frequency m^^) 



Hri)) = 



1 

1 + -H^{r] - T]^) 



During this period, S(f couples to the x-matter through the (^-dependence of the potential 
energy V{x, <f) ^ \m\{(p)x^ ^ |mJ((^o)x^ [l + iPx^V^] which yields 

VV - qV^^ix, p>o)Sp> = VV - '^PxQVix, Vo)5p> = 0. {A.22) 

Then the Fourier modes i^kiv) of Eq. (A. 4) satisfy the conformal-time evolution equation 

d'^Mv) + {k' - U{v))Mv) = 0, (^.23) 

U{V) ^ ^ - qaX'^ix, ^o) = ^- 2qP^a'V{x, <Po). {A.2A) 

a (X 

The first forms of equations (A.22) and (A. 24) are valid in both the infiationary period (where 
V^'ipiXi 'Po) = 4/?iF(x, (po)) and the x-dominated one. During inflation the effective potential 
reads U{r]) = 2(1 — 3/?i)?7~^, while during x-domination 

^(^) = _ 3^^)2 [1 - 3/5x + 3(1 - Px) cos[2m^(t - Q]] {v>V*)- {A.25) 
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The effective potential (A. 25) contains two distinct spectral features: a monotonic piece 
oc (?7 — 3?7*)~^ varying on the (averaged) Hubble time scale, and an oscillatory piece involving 
cos[2mj^(t — t*)]. Correspondingly, the spectrum of quantum fluctuations generated by solving 
the Schrodinger-like equation (A. 23) will have two peaks: one at A; ~ and one at A; ~ m^. 
[Depending upon the inflationary scenario considered, these peaks might be separated or 
might overlap.] One can compute the Bogolyubov coefficient of the first peak by matching 
at = ry* the exact solutions of Eq. (A. 23) when dropping the cosine term in Eq. (A. 25). As 



above, these are Hankel functions with v — = i^Q{j3i — 3/8) in the infiationary period and 



V = v.^ = iy/6(l3^ — 3/8) during the x-dominated one. (We assume here for definiteness that 
both are imaginary). This yields for the "Hubble-time scale" peak in the long wavelength 
limit k <^ 

1 I;/ I r / I \ ^ / 5^ \ ^ 

(yl.26) 



|2 
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Wx\ 


a- 








Ih 
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(4^x1) 



Note that this \Pk\ii is independent of k, and therefore different from the Rayleigh-Jeans-type 
spectrum (A. 14). Integrating (A.26) up to some effective cut-off Ajmax = leads to 

"^^^'^^ ^^-''"^ 

The value of k depends upon the details of the transition between inflation and x-domination. 
We expect that k,<1. Assuming k ~ 1, we get the total energy density (4.4) roughly inde- 
pendently of whether inflation is followed by a radiation-dominated era or a x-oscillatory 
one. 

Finally, the "oscillatory" peak k ~ can be estimated by applying Born perturbation 
theory to Eq. (A.23): namely, pa {2ik)~^ j U'^'"'{r])e~'^''^'^dr]. This yields {9 denoting the 
Heaviside step function) 

(o \ 2 rrS 3/2 

gj <(ix-^?^^0{k-m^l (A28) 



and, for instance, 

"r = ^(/'x - D^f . (A29) 

[These estimates are accurate when {(3-^ — l)H^ <S and should give the right order of 
magnitude when P^ — lr^l.] In order of magnitude the number density (A. 29) is comparable 
to (A. 27a). The corresponding energy density would differ from (A. 27b) by a factor ~ m-^/H^. 
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